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Summary 


The technique of implicit differentiation has been used in combination with 
linearized lifting-surface theory to derive analytical expressions for aerodynamic 
sensitivities (i.e., rates of change of lifting pressures with respect to general 
changes in aircraft geometry, including planform variations) for steady or oscil- 
lating planar or nonplanar lifting surfaces in subsonic, sonic, or supersonic flow. 
The geometric perturbation is defined in terms of a single variable, and the user 
need only provide simple expressions or similar means for defining the continuous or 
discontinuous global or local perturbation of interest. Example expressions are 
given for perturbations of the sweep, taper, and aspect ratio of a wing with trape- 
zoidal semispan planform. In addition to direct computational use, the analytical 
method presented here should provide benchmark criteria for assessing the accuracy of 
aerodynamic sensitivities obtained by approximate methods such as finite geometry 
perturbation and differencing. The present process appears to be readily adaptable 
to more general surface-panel methods. 
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Nomenclature 

A A 

coefficient of pressure-mode function f (£)g (h) 
(Eq. (10)) " m 

./?T7 


in pressure expansion 


b(y) streamwise semichord at spanwise station y 

b Q reference length (typically root semichord) 

AC complex amplitude of lifting-pressure coefficient 

P 

A 

f n (£) nth chordwise pressure-mode function 

A 

g m (n) mth spanwise pressure-mode function 

I j(kj),Kj(k 1 ) modified Bessel functions of first and second kind, respectively 
kernel function 

kernel functions for sidewash and upwash, respectively 
* b 0 w/U, reduced frequency 
- kr Q * kr/0, frequency parameter 
modified Struve function 

arc length of lifting surface, root to tip, in cross-stream direction 
free-stream Mach number 
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r All coordinates are normalized with respect to b Q . 
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x,y,z 


perturbation of lifting surface 

, 2 . 2\\/2 
- (x q + r ) 


Br, 


2 ^ 1/2 


- (y + z Z ) 
w o o' 

free-stream speed (in X direction) 

■ (MR - x Q )/6r 

complex amplitude of normalwash at point on lifting surface 
fixed streamwlse, lateral, and vertical coordinates, respectively 
Cartesian coordinates 


x o»Yo» z o x " 5 * y~ n » z_c » respectively 
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strearawise location of mean chord line 
streamwise location of leading edge 

A A 

coordinates on lifting surface transformed so that x - € « -1 at leading 
edge, x ■ 5 * +1 at trailing edge, and y ■ n ■ ±1 at tips 

* A - M 2 

lateral inclination (local dihedral) of lifting surface 
sweep angle of lifting-surface midchord line 
Cartesian coordinates 
- r/B 

circular frequency of oscillation 


Introduction 

Accurate and efficient computation of aerodynamic sensitivities (i.e., rates of 
change of surface pressures and weighted integrals of pressures with respect to 
changes in aircraft geometry) is of growing importance for aerodynamic shape optimi- 
zation as well as for multidisciplinary design synthesis (ref. 1). Such sensitivi- 
ties can, of course, be obtained approximately with any aerodynamic method by simply 
differencing the results of a solution for a base configuration and a solution for a 
geometrically perturbed configuration. The resulting sensitivity values, however, 
are dependent upon the size of the perturbation (step size) selected and are subject 
to cancellation errors if the step size is too small. Calculating sensitivities by 
analytical methods avoids such problems. 

Some capabilities to generate sensitivities and to solve the inverse 
(aerodynamic-design) problem have been developed from surface— panel methods (e.g«, 
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refs. 2 to 6), but these have been limited to steady state, and applications appear 
to have been generally limited to geometric perturbations of lifting surfaces and 
only in the direction normal to the chord plane. 

This paper presents expressions for aerodynamic sensitivities derived from 
subsonic, sonic, and supersonic unsteady, nonplanar lifting-surface theory (refs. 7 
and 8) for general perturbations, including planform variations. The geometric 
perturbation is defined in terms of a single variable. Moreover, without loss of 
generality, perturbations are performed by varying only the normalwash distribution 
and/or the streamwise planform coordinates. Performing the perturbations in this way 
accomplishes considerable simplification of the resulting expressions and avoids a 
third-order singularity along the line r * 0. Note that aspect ratio and even the 
spanwise location of a leading-edge or trailing-edge break point can be varied by 
perturbing only streamwise coordinates. The present formulation requires only that 
the user provide simple analytical expressions or similar means for defining the 
continuous or discontinuous global or local perturbations of interest (e.g., pertur- 
bation of sweep angle). In contrast, the method of reference 6 requires that the 
perturbation be made up of individual perturbations of each of the panel nodes - a 
computationally expensive operation. 

In addition to the steady-state limiting condition, the expressions derived here 
can be used, for example, to track flutter and other dynamic-response characteristics 
in multidisciplinary design processes. Such computations of aerodynamic sensitivi- 
ties, both steady and unsteady, should be useful in both conceptual and preliminary 
design. On the other hand, lifting-surface-theory calculations are themselves quite 
fast so that computation of sensitivities by differencing is not prohibitively expen- 
sive. Even so, the analytical expressions presented here provide convenient means 
for computation as well as accurate benchmark criteria for assessing the accuracy of 
differencing (e.g., as a function of perturbation size) so that, when needed, differ- 
encing can be used with better understanding and confidence. Hopefully, the present 
development will serve heuristic purposes as well. 


Theoretical Development 


Subsonic Flow 

The integral equation relating lifting-pressure and normalwash distribution on a 
nonplanar lifting surface (Fig. 1) oscillating in subsonic flow (Eq. (1) of ref. 7) 
can be written 



W N (x,y,z) 

_ 


j I AC b(n)K d£ dn 
-1 -1 P 


where 


K - Ky cos[y(y) - Y(n)l - Ky sin[y(y) - Y(n)1 


( 1 ) 


and Ky and K w are given, respectively, by 
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Ky(kx 0 , ky 0 , kz Q , M) - e ° 3 2 E/3y3z 
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and 
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E(kx , kr, M) - / * du 

u i J\ * u 2 


-ikjU 

VV • 1 7 'W - W> • ^ *--- *■ 
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(see Eqs. (4) to (6) of ref. 7). Carrying out the differentiation Indicated in 
Eqs. (2) and (3) leads to 
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p i (k i ) - KjCkj) + i y flj (kj ) - L^kj)] - i 


(7) 


Qj(k!> - K^kj) + k^Ckj) + i j [I^kj) - k 1 I Q (k 1 ) 


- L^kj) + k 1 L <> (k 1 )} - i 


( 8 ) 
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(see Eqs. (8) to (11) of ref. 7). The physical variables are related to the 
transformed variables in Eq. (1) by 


x - x L (y) + (x + l)b(y) - x m (y) + xb(y) 


5 ■ x L (n) + (5 + l)b(n) - x m (n) + £b(n) 


y “ b“ y 


t ^ 

n ” b“ n 
o 


( 9 ) 


The normalwash W N is defined by the shape and/or motion of the lifting surface. 

It is customary to represent the unknown pressure distribution AC by a linear 
combination of chosen shape functions. Thus p 


AC p a,n) 



L_-_S v-i, 


I b“ 1 (ri) l l a f (C)g(n) 
1 + 5 n m nm n m 


( 10 ) 


and the basic aerodynamic sensitivity is 
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f (Og(n) 

n m 


do 


where Q represents the geometric perturbation of the lifting surface and should not 

-1 * 

be confused with the function Oj(kj) in equations (5) and (8). Now — ^ is 

da 

purely geometrical; consequently, attention is focused on evaluating — g^- . Moreover, 

without loss of generality, the perturbation 0 will be performed by varying only 
the normalwash distribution % and/or the streamwise coordinates x and C in 
order to avoid having to deal with spanwise derivatives of K v and Ky which in- 
volve a third-order singularity at r ■ 0. Variations of x and £ may be defined 
in terms of changes in b and x^ or x m (Eq. (9)). Note that even a geometrical 
parameter such as the spanwise location of a leading-edge break point may be varied 
by changes in x^(y) only. 
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By use of Eq. (10), Eq. (1) may be expressed as 


H a F 
L L nm run 


b W„ 
o_ Q jj 
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where 


F 

nm 


1 1 J * ^ 

j I A - n 2 4 f n <5>* m (h){Ky cos[y(y) - y(n)l 


- Ky sin(Y(y) - Y(n)]} di dn 


( 12 ) 


* 


(13) 


Solution of Eq. (12) requires, of course, that the equation be written for at least 
as many points x,y,z on the lifting surface as there are values of a nm to be 
calculated. Differentiating Eq. (12) with respect to 0 gives 


l l 



3F 


nm 

3Q 
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(14) 


in which the values of F nm and a nn| are obtained from solution of Eq. (12) for the 
base configuration, and 3(W N /U)/3Q is geometrically determined and known a priori. 
Equation (14) shows that if the perturbation Q does not involve a planform change, 
then ^ n|B /3Q ” 0, and the required values of ^ a nn/^ ^ obtained fro® a direct 

solution of Eq. (12) with the perturbation normalwash 3(W N /H)/3Q input in place of 
%/U. 

Differentiating Eq. (13) with respect to Q gives 


3F 

nm 

3Q 


/' / JT77 E t (b. (hfi coslT(y) - y(i)1 

-1 V 1 + c m 1*0 


Note that 


3Ky 

3x 


sin[y(y) - Y(n)l 
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( 16 ) 
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is the only quantity in Eq. (15) that depends on the particular perturbation 

Q. 


Geometrical Example 

The example of a planar wing with trapezoidal seiuispan planform serves to 
illustrate the nature of 3x/3Q for perturbations in aspect ratio A, taper ratio 
A, and midchord sweep angle A. The semispan A is invariant. The first of 
Eqs. (9) gives 


x ■ 1 + y| (1 + i) tan A + x[l - (1 - A)y] 


Then with A and A invariant, 


3x _ 3x + A^ k 

■55 1 1a ■ tsn A 

With A and A invariant, 

3x _ 3x */ *. A . \ 

■55 * IX - n x + 2 tan A J 

With A and A invariant, 

3X _ 3X “A,. , , « , , , _ 2 ks 

Iq =TK " y 7 ° X)(1 + tan A) 

For use In Cq» (15) these expressions lead to 

3 x i \ 

o / A * v 1 + A ^ . 

1a“ ” 2 tan A 

3X A A A A A A ^ 

yy * yx - n£ + (y - n) J tan A 

- (y - n) j (1 + A ) ( l + tan 2 A) 
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and these may be combined into a single expression for use in Eq# (15). 


3x A A A A A A 

lo 2 " c o yx + c l (y ” n) + c 2 n5 


where for Q ■ A, 

c Q - 0, Cj ■ - 2 ~ tan A» c 2 “ 0 

for Q * X, 

A 

c o * C 1 " 7 tan A » c 2 “ _1 
for 0 ■ A, 

c Q * 0, Cj ■ j (1 + A)(l + tan^A), c 2 * 0 


Derivative of the Kernel Function 

The required expressions for 3Ky/3x 0 and SKy/dxQ in Eq. (15) are obtained by 
differentiating Eqs. (5) and (6), respectively. Thus 


3K V 

3x 

o 




(17) 


+ 



- lkK w 


(18) 


The singularities of K y and Ky (Eqs. (5) and (6)), and in particular the 
r singularity, remain in 3Ky/3x Q and 3Ky/dx Q (Eqs. (17) and (18)) only for 
unsteady flow. For both steady and unsteady flow, however, 3K w /3x {) contains a 
third-order singularity (R“ 3 ) at the point x Q * y Q * z© “ As indicated earlier, 
a third-order singularity along the line r Q - 0 has been avoided by not perturbing 
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in y direction. The R" 3 singularity may be dealt with in Eq. (15) by the 
classical analytical extraction technique. The remaining numerical quadrature, 
however, converges slowly. Results can be improved if second-derivative terras are 
also included in the extraction. 

Eqs. (17) and (18) contain, respectively, the factors y Q z 0 /R 5 and z 3 /R 5 
which also appear to indicate third-order singularity. However, if the normalwash 
control points (x,y,z) are not located at points of normalwash discontinuity (as, in- 
deed, they should not be!), then as y 0 goes to zero, z Q behaves as y^ (Rig. 1). 
Consequently, y 0 z o^ R ^ behaves as R”^, and z^/R 5 behaves as R“*. Thus these 
terras require no special attention in Eq. (15) beyond the aforementioned extraction 
of the third-order singularity. 

The use of Eqs. (15) to (18) permits Eq. (14) to be solved for the required val- 
ues of 3 a nm /3Q which, in turn, permit the evaluation of 3AC p /3Q from Eq. (11). 

If this procedure is employed for a number of different perturbations, only the sim- 
ple geometric expressions 3b _1 /3Q in Eq. (11) and 3 x q / 3Q in Eq. (15) need be 
reevaluated for each. 


The expressions for K v and Ky as well as for 3K v /3x 0 and <*Kw/3x 0 simpli- 
fy considerably for steady flow (k ■ 0) or for planar geometry (z Q - 0). Thus for 
steady flow, k * 0, and 


and 



so that 
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k*0 
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For planar geometry, z Q - 
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so that 
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Finally, for planar steady flow, k • z ■ 0, and 


and 


*W 


k»z -0 
o 




8lt w 



k-z *0 
o 



(26) 


(27) 


Alternative Formulation for Subsonic, Sonic, 
and Supersonic Speeds 

An alternative form of the kernel function has been used to define the kernel 
for sonic and supersonic speeds as well as for subsonic (ref, 8), Thus in the 
present notation 


K * e °(K 1 T 1 /p 2 + K 2 T 2 */p 4 ) (28) 

where 

T l - cost Y(y) - Y(n)l (29) 
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( 30 ) 


t 2 * ■ [z 0 cos y(y) “ y Q sin y(y)][z 0 cos y(n) - y Q sin y(n)l 
Kj - p(3l 0 /ap) (31) 


K 2 - p 3 (3/3p)[(l/p)31 0 /3p] 


and 


s 0 -iks 

i.- l 2 -j «■ 

s i J7T? 


(32) 


(33) 


Kj in Eq. (31) should not be confused with the modified Bessel function Kj(kj) in 
Eqs. (7) and (8). The limits of integration s^ and s 2 are determined by the time 
required for transmission of small disturbances from sending point (€,n,?) to receiv- 
ing point (x,y,z) with consideration of the limits imposed by Mach lines, if any. 

The derivative 3K/3x Q required for the generation of aerodynamic sensitivities 
may be obtained from Eq. (28). 
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p 2 + s 2 



(36) 
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The derivatives of sj and s 2 indicated in Eqs« (35) and (36) may be obtained from 
the expressions for Sj and s 2 that are pertinent to subsonic, sonic, or super- 
sonic flow (ref« 8). Thus for subsonic flow 



s 2 ■ °°, and all the derivatives are zero as in Eq. (38), 
For supersonic flow, 
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9s, 


_ o o 9 i 

>! - (x Q - MR)/B 2 “ (M p - x o )/(MR + x q ), - Mp/R, 


'N 


3s l\ 2 2, 2 8 S 1 h/ B 


S’ ^■( I - M r)/' 2 . 


3x 


■(if) ■ - m « 0 /r3 > ir(-5r) --2« 2 pV r4 > 


a 2 

3 I 3 S I 


x / a 2 N 

« , 9 , - -m -4 3ir -4 + i 

^oW 2 / r 3 \ r 2 > 


( 40 ) 


o 9s 2 / 9s ?\ 

s 2 - (x Q + MR)/B 2 , Tp Mp/R, {—) = M 


V/R 2 , «^(i+b 2 4' 

8p 2 R \ R 2 / 




IF - (' + M ir)/” 2 - it(-5f) - M "o /R ’- Mtt) • - 2mV V r4 > 


(41) 


9 ( 3 tl \ _ L 2 / ' 

9x o\3p 2 / R 3 V R 2 \ 


Substitution of the appropriate expressions from Eqs. (37) to (41) into Eqs. (35) 
and (36) and the latter two equations into Eq . (34) produces the kernel-function 

derivative needed for subsonic, sonic, or supersonic speeds. The kernel derivative 
expression is used, in turn, to evaluate 3F nm /90 from 


9F 1 1 i ^ ~a * R ~ Z a a 9x * A 

Jr—j f n < 5 V R) Ir 7 <r dC dn (42) 

-1-1 v 1 + £ o 

which corresponds to Eq. (15). 

Note that the weight factors (two radicals) in Eq • (42), as well as in 
Eqs. (10), (11), (13), and (15), are appropriate for subsonic speeds. Different 

factors and different pressure-mode functions f (£) and g (n) are used for sonic 

n m 

or supersonic speeds (see, e.g., ref. 9), and the domain of Integration is limited to 
the portions of lifting surface lying within the Mach forecone with vertex at the 
receiving point (x,y,z). 
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Concluding Remarks 


The technique of implicit differentiation has been used in combination with 
linearized lifting-surface theory to derive analytical expressions for aerodynamic 
sensitivities (i.e., rates of change of lifting pressures with respect to general 
changes in aircraft geometry, including planform variations) for steady or oscil- 
lating planar or nonplanar lifting surfaces in subsonic, sonic, or supersonic flow. 
The geometric perturbation is defined in terms of a single variable, and the user 
need only provide simple expressions or similar means for defining the continuous or 
discontinuous global or local perturbations of interest. Example expressions are 
given for perturbations of the sweep, taper, and aspect ratio of a wing with trape- 
zoidal semispan planform. In addition to direct computational use, the analytical 
method presented here should provide benchmark criteria for assessing the accuracy of 
aerodynamic sensitivities obtained by approximate methods such as finite geometry 
perturbation and differencing. The present process appears to be readily adaptable 
to more general surface-panel methods. 


Acknowledgement 

The author extends his sincere appreciation to Robert N. Desmarais for checking 
the derivation of equations and for valuable and stimulating discussions relative to 
numerical implementation. 




14 



References 


1. Sobieszczanski-Sobieski, Jaroslaw: The Case for Aerodynamic Sensitivity 

Analysis. In "Sensitivity Analysis in Engineering." NASA CP-2457, 1987. 

2. Bristow, D. R.; and Hawk, J. D.: Subsonic Panel Method for the Efficient 

Analysis of Multiple Geometry Perturbations. NASA CP-3528, 1982. 

3. Bristow, D. R.; Hawk, J. D.; and Thomas, J. L.: Subsonic 3-D Surface Panel 

Method for Rapid Analysis of Multiple Geometry Perturbations. AIAA 

Paper 82-0993, 1982. 

4. Hawk, J. D.; and Bristow, D. R.: Subsonic Surface Panel Method for Airframe 

Analysis and Wing Design. AIAA Paper 83-0341, 1983. 

5. Bristow, D. R.; and Hawk, J. D. : Subsonic Panel Method for Designing Wing 

Surfaces From Pressure Distribution. NASA CR-3713, 1983. 

6. Hawk, J. Dennis; and Bristow, Dean R.: Development of MCAERO Wing Design Panel 

Method With Interactive Graphics Module. NASA CR-3775, 1984. 

7. Yates, E. Carson, Jr.: A Kernel-Function Formulation for Nonplanar Lifting 

Surfaces Oscillating in Subsonic Flow. AIAA Journal, Vol. 4, No. 8, Aug. 1966, 
pp. 1486-1488. 

8. Harder, Robert L.; and Rodden, William P. : Kernel Function for Nonplanar 

Oscillating Surfaces in Supersonic Flow. Journal of Aircraft, Vol. 8, No. 8, 
Aug. 1971, pp. 677-679. 

9. Richardson, J. R.t A Method for Calculating the Lifting Forces on Wings 
(Unsteady Subsonic and Supersonic Lifting-Surface Theory). R. & M. No. 3157, 
British A.R.C., 1955. 


15 



IWVSA 

Natona* Aenonauhcs and 
Space Admrustraion 


1. Report No. 

NASA TM- 100502 


4. Title and Subtitle 


Report Documentation Page 


2. Government Accession No. 


3. Recipient's Catalog No. 


5. Report Date 


Aerodynamic Sensitivities From Subsonic, Sonic, and September 1987 

Supersonic Unsteady, Nonplanar Lifting-Surface Theory e. performing organization code 


7. Author(s) 


8. Performing Organization Report No. 


E. Carson Yates, Jr. 


10. Work Unit No. 


9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665-5225 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 

15. Supplementary Notes 


505-63-1 1 

11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Memorandum 


14. Sponsoring Agency Code 


16. Abstract 

The technique of implicit differentiation has been used in combination with 
linearized lifting-surface theory to derive analytical expressions for aerodynamic 
sensitivities (i.e., rates of change of lifting pressures with respect to general 
changes in aircraft geometry, including planform variations) for steady or oscil- 
lating planar or nonplanar lifting surfaces in subsonic, sonic, or supersonic flow. 
The geometric perturbation is defined in terms of a single variable, and the user 
need only provide simple expressions or similar means for defining the continuous or 
discontinuous global or local perturbation of interest. Example expressions are 
given for perturbations of the sweep, taper, and aspect ratio of a wing with trape- 
zoidal semi span planform. In addition to direct computational use, the analytical 
method presented here should provide benchmark criteria for assessing the accuracy 
of aerodynamic sensitivities obtained by approximate methods such as finite 
geometry perturbation and differencing. The present process appears to be readily 
adaptable to more general surface-panel methods. 


17. Key Words (Suggested by Author! s)) 

Aerodynamic sensitivities 
Unsteady aerodynamics 
Lifting-surface theory 


18. Distribution Statement 


Unclassified - Unlimited 
Subject Category - 02 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21 . No. of pages 

22. Price 

Unclassified 

Unclassified 

17 

A02 


NASA FORM 1626 OCT 86 

















